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Abstract 



^ l' We present approximations of varying degree of sophistication to the 

integral equations for the (gluon) structure functions of a hadron ("the 
partonic flux factor" ) in a model valid in the Leading Log Approximation 
with a running coupling constant. The results are all of the BFKL-type, 
i.e. a power in the Bjorken variable Xg^ with the parameter A deter- 
mined from the size uq of the "effective" running coupling a = 3as/n = 
ao/log(fc^) and varying depending upon the treatment of the transverse 
momentum pole. We also consider the implications for the transverse 
' momentum {k±) fluctuations along the emission chains and we obtain an 

exponential falloff in the relevant k = log(fc|^)-variable, i.e. an inverse 
power with the same parameter A. This is different from the 

BFKL-result for a fixed coupling, where the distributions are Gaussian in 
the K-variable with a width as in a Brownian motion determined by "the 
length" of the emission chains, i.e. log(l/a;s). The results are verified 
by a realistic Monte Carlo simulation and we provide a simple physics 
motivation for the change. 
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1 Introduction 



1.1 The Contents of this Note 

In an earlier paper, Q, we have presented the Linked Dipole Chain (LDC) 
model as a generalization of the well-known CCFM model (for Ciafaloni-Catani- 
Fiorani-Marchesini) , j^, to describe the hadronic structure functions in Deep 
Inelastic Scattering (DIS) events. We have further in [|| described a set of results 
from the model. In this note we would like to continue the investigations and 
in particular describe the solutions to the equations for the structure functions 
when a running coupling is introduced. We will find power-solutions in the 
Bj or ken- variable, i.e. the (gluon) structure function behaves as in the same 
way as for the BFKL-solutions Q for a fixed coupling. 

The power A is determined by the strength uq in the effective running cou- 
pling a = Sas/n — Q;o/ln(fc^) (In a pure Yang-Mills theory we have used the 
value ao — 12/11.) and the treatment of the transverse momentum pole, i.e. 
the small virtualities in our equations. Actually, it will turn out that the equa- 
tions become unstable in the sense that the value of the parameter A is sensitive 
to the approximate treatment of the running coupling in the soft region where 
In(fc^) is small. 

If we introduce a cutoff in the ln(fc^)-variable we obtain, however, inside a 
large region, very stable results for an isolated largest eigenvalue A. One major 
finding is that the LDC model in itself contains a suppression of the soft region 
close to the pole of the coupling constant, and thus contributes significantly to 
this stabilization (This stems from the averaging over the azimuthal angles). 
Our result is then for the value of the isolated largest power A that it is about 
A ~ (0.3 — 0.4)q;o within a realistic cut-off region. 

For the transverse momentum dependence of the structure function we will 
find the asymptotic solution (ln(fc5^))""/'^. Further, the transverse energy distri- 
bution along the chain will not, as for the BFKL case with a constant coupling, 
be Gaussian in log(/c^) with a width determined, as in a Brownian motion, 
from the "length" of the emission chains, i.e. log(l/xB) Instead it will be a 
(negative) exponential in log(/c^), i.e. it will behave as an inverse power in the 
transverse momentum kj^^^^^^ with the same A-value as above. (Note that the 
Rutherford parton scattering will behave as kj*. The extra factor kj^^ is due 
to the larger x^-values needed for higher k±.) This is also born out by a Monte 
Carlo simulation. 

It is a well-known property of diffusion equations with a force (in this case 
stemming from the running coupling which favours small fcj^-values) that there 
are on the short time scale (for small log(l/a;B)-values) an erratic behaviour, 
which can easily be mistaken for a stochastic Brownian motion, before the pro- 
cess reaches its long-time stable distribution. Our estimates of the scales in this 
case unfortunately indicate that the HERA range is too small for a noticeable 
change to the above-mentioned power behaviour in k±. 

We will end this note with a set of simple examples to show the way the 
model interpolates between the DGLAP |^ and the BFKL mechanisms, the 
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Figure 1: Fan diagram for a DIS event, qi denote the emitted ISB partons, ki 
denote the propagators and the dashed hnes are the FSB partons. 



reasons why and the mechanisms by which the stable distributions in transverse 
momentum emerge. 

1.2 The CCFM- and the LDC-Models 

To introduce the models we note that it is always possible to subdivide the 
radiation in the states into one part, the Initial State Bremsstrahlung (ISB) 
(denoted by the vectors qj in Fig |l|) and another part, the Final State Brems- 
strahlung (FSB) (the dashed lines in Fig |l|). The main requirement is that the 
FSB partons should be possible to emit in accordance with the QCD coherence 
conditions and with negligible recoils, if the ISB partons are already emitted. 
The observable emission weights are given by 



dw = J2dw^°\l)Ml) (1) 

i.e. by the sum over the chosen ISB sets (/) with the basic weig 
to emit these particular partons and a form factor A(/) stemming from the 
radiative corrections from the choice. To be precise each (/)-state contains a sum 
over all the exclusive states, containing these ISB partons and any other FSB 
partons. The FSB should have the property that the corresponding form factors 
Sud{F) are of the Sudakov type, meaning that, for each fixed (/)-state, the sum 
over all possible exclusive final states F becomes unity J^f du!i{F)Sud{F) = 1. 

If we ask for a particular value of the observables {Q^, xb), i-c. the (squared) 
"inverse wave-length characteristics of the probe" and the scaled energy-momen- 
tum, the Bjorken variable xb, then we should in Eq (Q) sum over the ISB sets, 
which end on such a configuration. 

The choice of the ISB set in the CCFM model is to order all emissions in ra- 
pidity (due to the relation between angle and rapidity this means angular order- 
ing, a well-known formulation of the QCD coherence conditions). Then CCFM 
chose as ISB the emissions with the property that there is no emission further 
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along the rapidity ordering with a larger light-cone energy momentum qj^. The 
ordering is done from the target-hadron side, i.e. the target is supposed to have 
a (large) positive light-cone energy-momentum P = {P^,P^ ~ 0,0^). In the 
probe-hadron cms the probe (with —q^ = Q^) has a large energy-momentum 
along the negative light-cone q = {—q+, q-,0±). This choice of the ISB is then 
consistent with QCD coherence but not with the fact that there should be sym- 
metry between the hadron and the probe directions, 0. 

Assuming that the recoils can be neglected, the ISB partons with momenta 
qj are on their mass shells. Energy-momentum conservation at every ISB vertex 
implies that the connector vectors ( "the propagators" ) kj , cf. Fig |l|, fulfill 

j 

kj=P-Y^ qj ie kj = kj^i - qj (2) 
1 

CCFM introduce the variables {zj,q±j) such that k-^-j — 2;jfc_|_(j_i) and k±j = 
k±(j-i)—q±j- In the Leading Log Approximation (LLA), relevant for the CCFM 
model calculations, the fractional variables should fulfill Zj <C 1, i.e. the splitting 
functions are approximated by P{z) ~ 1/z. While the emitted gluons along the 
fan diagram in Fig ^ are treated as mass less, the connector vectors kj are space 
like with —k"^^ k"^^-. This relation is satisfied provided 

kij > Zjqlj (3) 

while smaller fcj^-values are suppressed. Then the weight and the form factor in 
Eq (g) is in the CCFM model 

_ dz dq^, A / 7 N r.\ 

a ^Ane(.z,k^,q±_) (4) 

z ql 

with the effective coupling a — Sas/ir and the non-cikonal form factor 

A„e = exp(-a ln(l/z) ln{kl/zqlj) (5) 

The LDC model choice in Ref Q] is to restrict the ISB emissions in the CCFM 
model to those which also fulfill 

q± j > min ( fc^ J , A;^ _ 1 ) ) (6) 

Gluons which do not satisfy this constraint can be included in the FSB set 
causing only small recoils. If the ISB set is restricted in this way, and we sum 
over the corresponding weights and form factors it is in Ref jl] shown that we 
obtain the same weight as in Eq but with the non-eikonal factor exchanged 
for 1. This is obviously a major simplification in order to calculate the cross 
sections according to Eq (|l|). 

The second simplification comes with respect to the emission of the FSB 
radiation. The FSB gluons are emitted as dipole radiation with the "original" 
dipoles spanned between the color- adjacent ISB gluons {qj,qj+i) and with the 
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Figure 2: (a) The shaded area denotes the forbidden region for k±j where 
q± < min(k_Lj, kj_j_i). For 0.5 < /cxj/fc^j-i < 2, the (/)-intcgral is Umited to a 
restricted region, (b) The azimuthal average h as a function of hi(fc^^/fc^, -i))- 



largest allowed FSB transverse momentum determined by the propagator virtu- 
ality — fcj — k^j- This means that we may for the FSB radiation make use of the 
weh- known Lund Dipole Cascade model as it is implemented in the Monte 
Carlo simulation program ARIADNE, 

The transverse momentum restrictions for the FSB emissions, i.e. that the 
allowed FSB transverse momentum is limited by the propagator — fc| ~ fc^^ , oc- 
cur also in the CCFM model and is a major result from their very complicated 
calculations, |]. Physically it means that the two color adjacent gluon currents 
from the "dipole pair" (qj, qj+i) in the LDC model do not stem from the same 
space-time point. There is a distance, which due to Lorenz contraction is essen- 
tially transverse, b± ^ 1/ ^ —k'j. The constraint then stems from the difficulty 

to emit FSB radiation with a wave-length X ~\/q± < b, i.e. smaller than the 
"antenna size". In other words there should be a form factor, which if it is e.g. 
an inverse power in k± is exponentially falling in the relevant ln(fc^)-variable 
and therefore leads to negligible contributions to the LLA outside this region. 

After azimuthal (0) averaging for fixed |fcj^j_i| and \k±j \ over the transverse 
pole in = {k±j — fc^Q_i))^ we obtain with the constraint in Eq (H) (see Fig 
la) the LDC weights 

dwiLDC) = f^^f^e{q^,-min{k^„k^^,_,))) (7) 



a 



dzj dkj^ 
Zj max(fc^^., /c^|.^._-^-|) 



h{k\j/k\i^^_ 



The function h{a), which is obtained from the azimuthal integration, is given 
by (see Fig |b) 

, (l-f^^^tan(l±^y^)g(a-l/4)) 
/i(o) — ; < a < 1 
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h{a) ^ h{l/a) 



(8) 



As dzj/zj = dyj, with yj equal to the rapidity, we see that the expression in 
Eq. (n) is completely symmetric, and the chain could equally well be generated 
from the probe end towards the target. We also note that in the LLA in general 
ln(fc^^/A:^^^_-^^) is not close to and then h{a) ~ 1 (although we note that "at 

the pole" there is a dip h{l) = ^/3/2'k ~ 0.28). Thus in the LLA approximation 
(using K = \n{k'^) and £ — ln(l/z)) 



Therefore we obtain in the LDC model from Eqs (^ and (^ that the weight 
distributions are local (i.e. we may define a Markovian stochastic process). The 
locality property in particular implies that the process is simple to implement 
in a Monte Carlo generation program which simultaneously describes both the 
structure functions and the final state properties. Such a general program 
is now available and is linked to ARIADNE and JETSET 

The above-mentioned symmetry implies that we may use either the fractional 
variable z+j (from the target side) or z_ j (from the probe side) as variable in the 
Monte Carlo generation. In Refs and j|] we have discussed the possibility to 
extend the scenario outside the LLA. The result is that it is necessary to choose 
the direction, i.e. to use z+- or 2;_-generation, in accordance with the direction 
of increase in transverse momentum, i.e. in the virtuality, of the propagator. 
Then the use of the splitting functions in the proper light-cone fraction should 
contain many of the non-leading contributions to the structure function. The 
exact results for the radiative corrections are, however, not known in this case. 
Nevertheless the simplicity and symmetry of the result, i.e. the chain of linked 
dipoles (containing the FSB and then leading to hadronization in a similar way 
as in e^e^-annihilation events) with the restriction of no emission above the 
corresponding propagator virtuality, seems physically appealing and stable. 

In Fig 1^ the fan diagram emissions in Fig |l| are described in the well-known 
Lund dipole phase space (which for any dipole is a triangle in {y, k) with y the 
rapidity and with the height and the base-line size given by ln(MjjpQ;g)). The 
whole fan emission must occur in between the lines corresponding to ln(P_|_) 
for the target and ln{q-) for the probe (fixed values of the light-cone energy 
momenta In(fc-i-) — l/21n(fc^) ± y correspond to straight lines in the triangle). 
The emitted gluons are described as extended triangular folds, starting at the 
{y, K)-value of the on-shell gluon. The dipoles are spanned between the "tips" of 
these adjacent gluon folds, i.e. they again correspond to triangles (the two sides 
of the triangular folds correspond to the colour and anti-colour of the gluon, 
which are connected each to the two dipoles around the gluon corner). The 
emission region of each such dipole is limited by the maximum virtuality, which 
corresponds to the propagator ln(— fc^) ~ In(fc^). 

A propagator in this plot is specified by a horizontal line indicating its re- 
value (transverse momentum) and with the "starting point" (right) given by 




if Kj > Kj-l 

1 ) otherwise 



(9) 
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Figure 3: The Dipole phase space in DIS. The shaded area is the phase space 
for FSB. 

the rapidity value y- = \n{—k±/k-) and the "end point" (left) given by y+ = 
ln(fc_|_/fc^). This means that the difference 5y between the "endpoint" of one 
propagator and the "starting point" of the next is 

< <5y = ln(fci,„„J - ln(fci^„) (10) 

For completeness we mention that inside the LDC model it is possible, Ref 
[0, to consider different sub- processes in DIS. Besides the "ordinary" parton 
model description with exceeding all the other momentum transfers in the 
chains, there is the boson-gluon fusion process, where the the last splitting ex- 
ceeds and all the other ones. Finally in a Rutherford parton scattering pro- 
cess the largest momentum transfer, fc^maa;' ^^^^ occur somewhere in the middle 
of the chain. Then we may consider the result of the dipole chain generation 
as the convolution of the target and the probe structure functions (each devel- 
oped to the virtuality —kf-^^ax with the Bjorken variables x±, respectively) and 
multiplied by the well-known Rutherford transverse momentum pole (A:^a£c)~^- 
This also implies that such a Rutherford scattering k±^ must always be larger 
than all bremsstrahlung fcx in the event. This constraint produces an effective 
cut-off for small k±_ Rutherford scattering. 

2 The LDC Evolution Equations 
2.1 Integral Equations 

We study in this paper only purely gluonic chains, and we let the gluonic struc- 
ture function F{x,Q'^) denote the gluon density in ln(l/a;). (Thus F/x is the 
density in x.) In this section we will treat the LDC evolution equation for the 
"non- integrated" structure function, !F{x,k'j_^) defined by 
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Here is the transverse momentum of the last link in the chain. We note that 
the contribution from chains with k\ > is suppressed by the factor Q^/k 



2 

J.' 

2 /n2 



which makes it sub-leading. We also note the modified argument, x — > x-k^/Q 
obtained for this term within the LDC model. Although less important this 
modification further suppresses this term. 

From the weights in Eq (|^) we find a recursive relation for and in ||l| we 
derived the following evolution equation: 



-{£,k) = / dK'a{K)h{K- k')T{£,k') 



di 

+ / dK'a{K')h{K' — k) exp[~{K' ~ k)]T{£ + K ~ k' , k') (12) 



We have here used the variables £ — ln(l/a;) and k — ln(A;^), introduced in the 
previous section. The function h defined in Eq (^) originates from the azimuthal 
integration in Eq (^. From now on we express h in the logarithmic variables 
K and k', and note that it depends only on the difference |k — (Note that 
in the LLA we can put h = I. The argument of J- in the second integral is 
explained in Eq (p^)). 

The first integral in Eq ( p^ ) corresponds to chains, where the last step is 
upwards in transverse momentum from k' to k, and the second integral, where 
k' is larger than k, to chains with a final step downwards. For the scale in the 
running coupling we have taken the largest virtuality in each vertex. Thus in 
the first integral we have q;(k), while in the second term we have «(«;'). 

With a running as the result will (as we will see in the following) be rather 
sensitive to the behaviour in the soft region for small k±_ . Naturally a perturba- 
tive calculation cannot be trusted when as becomes very large. Lacking a good 
understanding of the soft region, some kind of cut-off is necessary, and in this 
paper we have studied the results obtained assuming a lower cut-off, k±cut or 



Kcut = l^(^icut/^ )' (|llD and (12), while keeping the form a — ao/n 



for the coupling constant. (Another possibility would be to study the result 
assuming that as saturates for fcj^ -values below some cut.) 

After a Mellin transform in xb, which corresponds to a Laplace transform 
in £ with 



fx{K) ^ J d£exp{-\£)T{£,K) (13) 
we obtain (we will in general neglect the index A on / from now on) 



7 



A/(k) = J dK'a{n)h{n- k')}{k') 

+ I dKa{K)h{K - K)eyiY>[-{\ + l){K - K)\f{K) (14) 



It is sometimes useful (of section U) to introduce a symmetric non-integrated 
structure function Tg by the definition 

k) = T{1, k) exp(-K/2) (15) 

In that way the exponential suppression factor exp— (k' — n) in the second 
integral of Eq (^2|), which is "paid" only "for going down" in k, is changed so 
that "we pay the same amount for going up as for going down" in k. 

2.2 The Case of a Fixed Coupling a 

The case of a fixed coupling has been treated in Ref iQ (and also at many 
other places) and we will only briefly describe (some of) the results for future 
reference. We perform a Laplace transform with respect to k, and define / by 

/(7) = J ^''^ exp(-7K)/(K) (16) 

As /i is a function of k — k', we actually work with a convolution integral, and 
pole-singularities are obtained in / at the 7-values (the "anomalous dimensions" 
for given A) which satisfy 

A = a[g(7)+g(A + l-7)] with 5(7)= / dA exp(-7A)/i(A) (17) 

Jo 

The two (/-contributions stem from the two integrals in Eq (p^, respectively, 
and h{ A) is here interpreted as a function of the logarithmic variable A = k — k'. 

There are different levels of approximations possible to apply in the treat- 
ment of the Eqs (12), ([1^ ) and (p7|). The simplest one is to make use of the 



approximation in Eq (W and to replace h hy 1, which gives g « I/7. The 
function h deviates from 1 around the singular point ln(K^/K|_]^) = 0, and in 
an improved treatment we can describe this deviation by a (5-subtraction (cf ref 
0): 

h~h=[l- AS{k - k')] (18) 
Inserted into Eq (|l^ this gives 

3(7) a [I/7 - A/2] ^ \c^a[- + A] (19) 

7 A -I- 1 — 7 
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If wc fix A by the condition /(/i(A) - l)dA = J (/i(A) - l)dA we find A ^ 
0.65. In g also other moments of h are important, and in ref ||^ it was found 
that for A « 0.8 the same A-value was obtained from the original h and the 
approximation h. 

We note that the BFKL result corresponds to a neglect of the A-dependence 
in the second g-teTui in Eq (^7|). Although this is a sub-leading term, it is 
numerically not small, as discussed in Ref 0|. (There is also in BFKL a different 
method to regularize the singularity at ln(Kj/Kj_i) ~ 0.) 

In ||l[ we also pointed out another important non- leading effect. For consis- 
tency it is necessary that "the steps" in ln(l/z) really are large, but it turns out 
that there are essential contributions to the BFKL-integrals from eg z > 1/2. 
We return to this problem in section 

The xs-dependence is then determined by the inverse Mellin (or Laplace) 
transform. The leading behaviour stems from a pinch singularity at A = As, 
when the two symmetrical solutions to Eq ( |l7| ) coincide, i.e. when 7 = (l-|-As)/2 
(or 7 = 1/2 if the A-dependence in g{X -|- 1 — 7) is neglected.) From Eq ( p^ 
we find e.g. for a = 0.2 and A = 0.8 the result As « 0.41 (neglecting the 
A-dependence in g, we would instead obtain As « 0.64), with the corresponding 
Iow-xb behaviour ~a;^^° . 

When the inverse Laplace(Mellin)-transform is performed one also obtains 
the dependence on k± (or rather on k = In(fc^)) by a saddle-point approximation 
around the pinch singularity. The result (besides essentially trivial kinematic 
factors) is a Gaussian in k with a width ~-\/ln(l/a;), a result which can be 
interpreted as a kind of random walk in K-space Q . 

2.3 The Case of a Running Coupling a 

For a running coupling, chains with larger k- values are disfavoured by the small- 
ness of Qfs . This will imply that < k > does not increase beyond any bound as 
in the case with a fixed coupling. Instead the K-distribution saturates, and for 
small a;, J^{x, k'^) approaches asymptotically a factorizing form 

T{x,k) ^ x-^f{K). (20) 

The result appears to be similar to the type of random walk an atmo- 
spheric molecule follows in the earth's gravitational field. In this case, when 
the molecules are constantly pulled downwards, an equilibrium is obtained for 
an exponentially decreasing density (see also the discussion in Appendix ^). As 
discussed in section |3.3| , also in our model with a running as the asymptotic 
fc_L-distribution is an exponential in k = In fc^ . Inserting the asymptotic form 



in Eq (|20| ) into Eq (12) implies that /(k) must satisfy Eq ( [14[ ) above, now with 
a{K) = ao/K, where ao = 12/11 if we study a purely gluonic situation with no 
quarks. The solution to this equation will be studied in the following section. 
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3 The Solution for Running Coupling 



We have not been able to find an exact analytic solution to the integral equation 
(p^. Instead we have solved it numerically, using Chebyshev polynomials to 
convert it to a matrix equation. This method was also used in ref and 
is discussed further in appendix The ^-dependence for small a:-values is 
determined by the largest eigenvalue A to this matrix equation. 

In order to better understand the properties of the model we have also 
studied a set of approximations, which all agree with the original Eq (^4|) to 
leading logarithmic accuracy: 

a) by using k instead of k' for the argument of a in the second integral of Eq 



(14) 



b) by using the approximation ft, « (1 — AS{k — k')) in Eq ([iq). A special 
case (called b') with A — corresponds to ft 1. 



c) with both of the approximations a) and b) (cf section 3.1). (Case c') 
corresponds to both approximations a) and b').) 

The cases b') and c) can be studied analytically. Case b') is discussed in 



appendix M, and the simpler case c) will be discussed in section 3.1 



3.1 Qualitative properties 

In order to understand the qualitative features of Eq (|lj) we start by investigat- 
ing the approximation called c) above, which is obtained when a{K') is replaced 
by a{K) in the last integral of Eq dlj), and the variation of h is approximated 
by a (5-function according to Eq (|l^). We note that both these approximations 
are of non- leading order, and the solution has properties similar to the solution 
of the original Eq (|lj). With these approximations the integral equation in Eq 
( p^ can be transformed by straightforward differentiations into the second-order 
differential equation 

u^^[{X + l)u-2]^^ + {X + l){l-^,)f (21) 

where u ^ k + fiA and fi ~ uq/X. This is equivalent to (the radial part of) 
the Schroedinger equation for the hydrogen atom with the angular momentum 
variable £ = 0. To see that, we introduce the function ijj defined by 

/ = exp[(A + l)/2)u]i:iu) (22) 

We obtain by straightforward calculations 

I d^u^/j) MA + 1)^ _ (A + 1)^V> 

2 udu^ 2u 8 ^ ' 

Thus the "mass"-value is unity and the (squared) "charge" equal to /i(A+ l)/2. 
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For large values of n (large values of u) the second term in Eq ( |2^ ) can be 
neglected, and the solution behaves as utp ~ exp[±(H-A)K/2]. Here the physical 
solution must correspond to the minus sign. From Eq (|2^) this implies that / 
varies as a power of k for large K-values. From the behaviour of eq. (^l|) for 
large u, ie large k, we find 

/ cx K^^^; K large (24) 

The value of A is determined by the boundary condition for small n, which 
can be derived from the integral equation ( [T^ ) . (With the present approximation 
we find /'(ftcut) = [1 + A — l/(Kc„t + ^A)]f{Kcut)-) This implies that we have 
solutions for a discrete set of eigenvalues A. The asymptotic behaviour of the 
solution is determined by the largest of these eigenvalues, which corresponds to 
a solution -0 with no zeros between Kent and oo. 

We note that since A is fixed by the boundary condition, the result is sensitive 
to the value of Kcut, i e it is sensitive to the soft region. When Kcut increases, 
the value of A decreases continuously. Note, however, that if we include the 
correction term in h then this sensitivity to Kcut is strongly reduced. As k = Kcut 
implies Ucut = Kcut + we keep away from the singular point u = in Eq 
(p3|), even when Kcut approaches 0. This stabilizing effect of h will be further 
discussed in the following subsection. 

For a fixed coupling we can from Eq ( p^ derive a differential equation similar 
to Eq (pl|). This equation does not, however, contain the factors u in front of 
(Pf/du and df/du. Therefore this equation has two exponentially growing 
solutions, which implies that in this case the dependence of upon x and k 
does not factorize. 

Summarizing wc have found that the asymptotic behaviour of is given by 
the form 

(xx-^ ■ K^-'^ = x-^ ■ (Inkl)^-^; ln(l/a;)>K>l (25) 

where the power A is an eigenvalue determined by the boundary conditions. It 
turns out that these qualitative features do not rely on the approximations used 
in this subsection, but are also relevant for the solution to the original equation 
®- 

3.2 The x-dependence 

We will now study the numerical solution to Eq (^4|). As the second integral is 
sub-leading, we see that for a = ao/re the eigenvalue A will to first approximation 
be proportional to uq. Since we are here studying only purely gluonic chains, 
and are not including quark links, it may be most consistent to use a value 
of ao which corresponds to a pure Yang-Mills theory. From the relation ao = 
36/(33 — 2nf) we find for a pure Yang-Mills ao = 12/11 as compared to e.g. 
ao — 4/3 for = 3. 

As we have not studied the influence from the quarks we present results for 
the more stable quantity A/qq. Thus in Fig (|^) we show how A/ao depends on 
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Figure 4: The dependence of X/ao on Kcut- The solid line is from cq. ( p^ 
without approximations. The other lines are for the approximations a, b' and c' 
from down and up. As the calculations do not include quarks we have presented 
results for the more stable ratio X/ao, where ag = 12/11 for Uf — and 
ao = 4/3 for Uf — 3. 



the value of Kcut in the range 0.01 < Kcut < 2 (To be exact the results presented 
correspond to ao = 12/11 and for ao = 4/3 they are about 4% smaller). Results 
are presented for the original equation (|lj) and also for the approximations a), 
b') and c') defined above. We see that the result is sensitive to both the soft cut- 
off, Kcut, and to the non- leading modifications in the different approximations. 
We note, however, that the factor h in the LDC model provides a very strong 
stabilizing effect. Thus for the original non-approximated equation, the result is 
fairly stable in a range 0.5 < Kcut < 2.0, which ought to contain realistic values 
for Kcut- In this range we find 0.3 ao < A < 0.5 ao. 

In the LLA the energy fractions Zi are assumed to be small compared to 1, 
and in ref we pointed out that quantitatively significant contributions arise 
from large z- values. We here show that this is true also with a running coupling. 
For the splitting function, we have in the weights in Eq (Q) used P{z) oc 1/z, 
which is correct within the leading log approximation. The full splitting function 
without regularization for z — I reads 

P(z) = 2A^,|i-^ + ^ + z(l-z)| (26) 

If inserted in an evolution equation, the pole at z = 1 has to be regularized. 
Virtual corrections imply that this pole does not contribute to the increase of the 
cross section (ie the structure functions) for increasing Q^. One way to interpret 
this is to note that in the gluon splitting process the "old" gluon disappears and 
is replaced by two "new" gluons. Thus only one new extra gluon is produced, 
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Figure 5: The dependence of A/ao on Kcut- The sohd hnc is for the default 
sphtting function. The dashed hne is for P{z > 0.5) = 0. As in Fig (Q) we 
present results for the more stable ratio A/ag. 



and values of z close to the pole 1/(1 — z) corresponds to the original gluon 
losing some of its energy. Thus this contribution does not increase the cross 
section, but only "shifts it" to slightly smaller x-values. 

In the LLA the energy fraction z, which corresponds to the "new" gluon, is 
assumed to be small. If we want to reduce recoils as much as possible, we can 
say that the "new" gluon by definition is the one with least energy, which must 
imply that z < 0.5 (c.f. refj . Therefore it might be sensible to replace the factor 
1/z in the weight in Eq (0) by {l/2Nc)P{z)0{O.5 - z). To get some estimate 
of this effect we note that P{z)/2Nc is smaller than 1/z in most of the region 
< z < 0.5. (Actually P{z)/2Nc~-l/z is positive only for z > 0.43, and also here 
it is fairly small.) Therefore to estimate qualitatively the consequence of this 
non- leading effect, we have studied the changes obtained from the replacement 
1/z — » l/z • 9(0.5 — z). It is straightforward to show that this change will give 
an extra factor 2^ to the left hand side of Eq (^4|), which means that A/ will 
be replaced by 2^ A/. Apparently, this leads to smaller values of A with an 
approximate factor 2^^ « 1/V2. In Fig ||, we compare X{Kcut) for the two 
alternative splitting functions. This non-leading correction is clearly significant, 
reducing A by 20-25%. 

We have also checked that the largest eigenvalue is well separated from the 
smaller ones (see Fig ^) . Thus for Kcut = 1 we find for the largest and the second 
largest eigenvalue the values 0.4 uq and 0.2 ag. Thus for ln(l/a;) > 5 the non- 
leading contribution should be supressed by a factor e^^ (and for ln(l/x) > 10 
by a factor e^^). 
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3.3 Transverse Momentum Dependence 

In Fig we show the sohition to the original equation together with ap- 
proximation b'), which corresponds to the replacement h 1. It is easy to 
demonstrate that in all cases /(k) and T have asymptotically the same pow- 



erlike behaviour as the simpler approximation studied in section 3.1, Eqs (^4|) 
and @. 

We note, however, that the parameter A takes on different values in the 
different approximations. If we calculate the first correction term to / in an 
expansion in 1/k 

f{K) (X K^-^l + bn-^ +0{k-^)) (27) 

we obtain 

1 



^ _ ao{ao - A) 
A 



l + A 



-A 



where A = 



h{A))dA « 0.65 (28) 



For A sa 0.4, b is very small, and from Fig ^ we see that the single power in 
Eq ( p^ ) is a surprisingly good approximation in the whole K-range. (For the 
approximation b') with h —^ 1, we have A — 0, and the correction term becomes 
much larger, as is also seen in Fig|^.) 

We note that the powerlike dependence on k for / and J- implies a similar 
power dependence also for the original "integrated" structure function F. For 
small X we have from Eqs (pO) and (El]) 



F = x-^F{lnQ^) 
F= f{K)dn + / e^i+^^^^'^e^ (29) 
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The result is presented in Fig which also shows the separate contributions in 
Eq (p9|). These contributions correspond to chains with fc^ < and > 
respectively, where k± is the transverse momentum of the last link in the chain. 
Inserting the asymptotic expression in Eq ( p7|) we find for large 



Fcx x-^(lng2)"«/^; ln(l/2;) >ln(Q2) > 1 (30) 

For large the first integral in eq (|2^) dominates, and for this term the single 
power in Eq ( |30| ) is a good approximation. The second integral is suppressed by 
the factor x{i+x) ' which for InQ^ ~ 10 is about 20% , but becomes relatively 
more important for smaller Q^. Thus, if we use the power expansion of / in Eq 
(p7|), we get for 6 « the better estimate 

F cx (InQ^) 



(31) 



To better illustrate the A;_L-dependence of the chain, we have studied a very 
long chain in order to find out where this chain passes the central line y = 0. 
It is straight forward to show that the density p (in the (y,K)-plane) of all the 
right hand endpoints of the horizontal lines in the diagram in Fig |3[ is given by 
the expression 

P = •exp(-K) (32) 

where £1 = InW — k/2 + y, £2 = InM^ — k/2 — y and W equals the total cms 
energy of the chain. The non- integrated function T(£2^ k) describes the weight 
for the chains coming from left in Fig ^ to the point {y, n), while the normal 
structure function F{£i, k) adds up the contributions from links with all possible 
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Figure 8: The In Q^-dependence of the structure function. The dashed line is the 
contribution from the first term in Eq. (29). The dotted hue is the contribution 
from the two first terms in the power expansion in Eq (^fj). 



values of k', both those larger than k and those smaller than k. Inserting the 
asymptotic expressions in Eqs (^5|) and (^o|) we find the following result (the 
factor exp(— Ak) comes from the ^-dependences of £i and £2) 



K^-ie^(i+-^)"; and K large (33) 

In Fig ^ we present the results from a Monte Carlo simulation for different 
energies W of the quantity p in Eq (^^. We note that the distribution has the 
expected approximate exponential fall off exp[— (1 + A)k], in accordance with 
Eq (^3|) and this corresponds to a (negative) power dependence in k± for the 
density p. Further the curve has the same slope independent of the total energy 
W, consistent with an asymptotically factorizing dependence on x and k±, i.e. 
there is no sign of the expected BFKL diffusion in k. The dashed line in Fig|^ 
shows the asymptotic form in Eq (|3^). 

For comparison we also show in Fig ^ the result for a constant coupling 
a = 0.2. Here we find a Gaussian distribution with a width increasing with the 
energy (~ -y/ln 1/x), as expected from a random walk. 



4 Dominating Paths in Transverse Momentum 

The result in Eq ( |30| ) is particularly interesting if we compare to a result obtained 
in 1^, where we study the mean paths of the initial state cascade in the {i, k)- 
plane. We find that due to the running coupling the region covered by the 
major contributions (nowadays known as the "Bartel cigar") corresponds to 
one part, length £1, with small k- values and with a BFKL-like contribution, 
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~exp(A£i) and a second part, length ^2 = ^ — ^i, corresponding to increasing n- 
values up to the final kq = ln((5^). This second part is a DGLAP contribution, 
~exp(2-y ao^2XQ), where XQ = 1ii('*q)- Using a saddle-point approximation 
we demonstrated that the combination of the two parts corresponds for > 
"oXq/A to 



which is identical to our result from Eq (pO|). 

For Xi < aoxq/X the DGLAP chains dominate, because in this case the 
saddle-point is not within the the allowed integration region. We note that 
this result implies that the BFKL behaviour is relevant only for rather small 
a;- values, x - 10"^ for - 100 GeV^ 

We will end by showing that for small x-values the very construction of the 
LDC chains leads to a simple result of the BFKL kind. In the LDC model it is 
possible to go "up" and "down" in k along the chains, and thus the total chain 
can be divided into "cells", where each cell contains a set of up-steps followed 
by another set of down-steps. Let us study one such cell, which starts in the 
point at (2/1, Ki), goes up to a maximum point (y', k'), and then goes down to 
the endpoint {y2,K2)- We note that this second part corresponds to up-steps 
from the probe side, and therefore both pieces correspond to "DGLAP-chains" . 

To describe the path we generalize the structure function to a function 
J^{£, K2, 1^1), which describes the sum of all paths, which start at a virtuality 
given by ki (which may now be different from Kcut) and end at K2 ■ It is also con- 
venient to include the exponential factor in Eq { p2[ ) into the structure function, 
and define a "symmetrized structure function" J-g by the relation 



If expressed in the rapidity separation between the two endpoints, Sy — £ + 
{k2 — Ki)/2, this function is fully symmetric with respect to the direction of 
the chain. The result for a long chain is also the direct product of the separate 
pieces. 

Inserting the DGLAP-result for the two sections in the cell, and integrating 
over possible intermediate points (y', k') we find 



We start by considering the simpler case with a fixed coupling a, for which 
we have 



(34) 



(35) 




(36) 



^DGLAP(£, k', Ki 



k' - Kl)) 



(37) 
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Inserted into Eq this gives 

Ts,ceu ~ J dy'dn' exp[2Va4(K'-Ki) + 2 V^I^I^T"^^ -k+ K2/2 + ki/2] (38) 

where €1 and £2 are given by Eq (^6|). This integral can be solved by stationary- 
phase methods, i.e. we look for the saddle-point in the integrand exponent. 
After a little algebra we find the surprisingly simple result that 

Ts,ceii ^ exp[{R- l)Sy/2] with i? = Vl + 8a and Sy ^ 1/2 - yi (39) 

Thus the contribution from one cell depends only upon the rapidity difference 
between the starting point and the endpoint. This also means that this result 
can be easily generalized to any number of cells. 

A closer examination tells us, however, that the maximum is only obtained 
within the integration region if Sy is large compared to 6k, = K2 — ki. There 
will be a dividing line with 

Sk - ^^Sy (40) 

(where R is defined in Eq (|3^)) with the property that for smaller k- values 
there is a saddle-point but for larger k the main contribution is a single DGLAP- 
motion always directed upwards in k. For k below the line the result is obviously 
of the BFKL kind, i.e. there is (besides the symmetrical K-dependence) an 
effective x~^^ behaviour, but this time with = {R — l)/2 ~ 0.3 for our 
"conventional" value ofa~3/7rx0.2. 

For a running coupling, cells which go very high up in k are suppressed. 
Therefore the cells will have limited sizes, both in k and in y, and for very small 
X the number of cells will be proportional to the total rapidity range, ln(l/a::). 
The result becomes a BFKL-like power-dependence on 1/x but a saturating 
distribution in k, which is described by the limited K-distribution in a single 
cell, a behaviour which is in agreement with our results in section |^. The 
dividing line in Eq ( ^ ) will be relevant in the region not too far from the target 
end, before the K-distribution saturates for very large rapidity separations. 



5 Conclusions 

We have shown that if wc introduce a running coupling a ~ ao/n into the 
integral equations for the gluon structure function in the LDC model then 

• for the non-integrated structure function J-{£ = ln(l/a;),K = ln(A:^)) we 
obtain a factorizing BFKL-like behaviour 

T{i,K)^eM^me)hA^^) (41) 

where the largest eigenvalue Am is isolated. It depends upon the soft cutoff 
Kc but takes on stable values A„i ~ 0.4 ao for 0.5 < Kc <2. Further a very 
good approximation for valid for k > Kc is 

/a„ ~ (42) 
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• For the (integrated) structure function F{£ = \n{l/x), kq = \n{Q'^)) we 
then obtain, again as a good approximation for large ^, {Xm£ > Q^o ln(KQ)/A,„) 

F{£, Kq) ^ const ■ exp{Xmi + ao In^Kg)/ A^) (43) 

which, as it is demonstrated in corresponds to an interpolation be- 
tween the BFKL and DGLAP mechanisms. For ao ln(KQ)/A„i > Xm^ the 
DGLAP mechanism dominates and there is evidently according to Eq ( ^3|) 
a smooth turnover to DGLAP when aoln(KQ)/A„i = Xm£. 

• The expected BFKL-diffusion in k (with £ as the "time" -variable) is not 
noticeable in the model. For large values of £ we obtain in the center of 
phase space a saturating ^-distribution (valid at least for k > 2) 

pdK ~ dKK^ao/A^-l ^ ^ ± - . (ln(fc^))2"o/A,„ (44) 

which may be characterized as a kind of "rcnormalizcd Rutherford cross 
section" . 



A Numerical Method for Solving the Integral 
Equation 

We will here briefly describe the method of using Chebyshev polynomials for 
solving the integral equation (p^). 

The Chebyshev polynomials are defined by 

Tn{x) — cosn(arccos x), — 1 < a; < 1 
and obey the orthogonality relation 

f.;-t'i„ ,45) 

J-i Vi-- X [ 7r,m = n = 

The eigenfunction and the kernel of the integral equation (times a factor 
\/l — x'^, because of the weight factor in the orthogonality relation) are ex- 
panded in Chebyshev polynomials, up to a certain degree N, thus converting 
the integral equation to an approximately equivalent matrix equation. Usually, 
we are interested in the largest eigenvalue and the corresponding eigenfunction. 

For practical reasons, we must approximate the upper limit in the integral 
equation with a large number which we have set to K^ax = 100- No significant 
dependence on the value of K,nax > 50 has been noted. 

Normally, the error in the largest eigenvalue A is less than 1% for the number 
of polynomials iV > 50 but we have used up to 320 polynomials to get a good 
description of the eigenfunctions for large k- values. 

In this work, we have used computer code from Numerical Recipes pM- 
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B Analytic Solution to the Approximated Inte- 
gral Equation 

In this appendix we wiU present an analytic solution to the integral equation: 

A/(k) = a r fiK')dK' + e^^+i)" T'"""' a/(K')e~^^+'^"'dK', (46) 
with a running coupling a{hi) — ao/n. This corresponds to approximation b') 



in section 3.2 where h is replaced by 1. 

For finite Umim i^max the integral equation is solved by rewriting it to a 
standard differential equation with solutions which are expressed in confluent 
hyper geometric functions. Using the boundary condition at Kmax and taking 
the limit K^ax ^ oo, we find that solutions to the integral equation have the 
asymptotic (k — > oo) behaviour /(k) ~ k~~^. Using also the boundary con- 
dition at Kmin, the largest eigenvalue can be found as a function of Kmin- No 
simple relation for \{Kmin) is found so it is evaluated numerically. The results 
are in agreement with the numerical solution of the integral equation. 

We start by changing variables to m = (A + 1)k to get rid of A from the right 
hand side of the integral equation ( ^ ) which becomes 

— /(«) = - / f{u')du' + e" / ^—f{u')du'. (47) 
This can be rewritten to the differential equation 



A 

\ w a / 

With 

yiu) = — ^— => ./(k) = /(m) = —y'iu), 

it can be rewritten to the following differential equation for y{u) 

1 , 1\ 
uj A 

Two linearly independent solutions are given by the real and imaginary part 
of the function (see e.g. Ref |l4[) 



y"-y'+{i2 + -.]7y-o. m 



M(u) = — ^ ^u^+^^M i-^v^ + w,l + i2v; u 



where = y ^ — j and M is the standard notation of the confluent hyper 
geometric function: 



M{a,b; u) = Y, 



a(a + 1) • • • (a + n - 1) m" 
b{b + l)---{b + n-l)~^ 
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The real and imaginary parts of M have the asymptotic behaviour 



yi{u) = Re M{u) 
1)2 {u) = Im M (u) 



-20. ^ 

u ^ e 

u >■ 



(50) 



To find the boundary conditions from the integral equation we insert the 
differential operator in the left hand side of eq. (^8|) into the integral equation. 
The boundary conditions are 



/' - / 



= 



1 



(51) 



From the second boundary condition, we see that the power solution (j/2 in 
Eqs. (pO|)) becomes more and more important as Umax oo- This means that 
solutions to the integral equation have the asymptotic behaviour /(k) ~ . 

Having singled out ?/2, we use the lower boundary condition, which simply 
means y(umin) — 0, to find a relation between the largest eigenvalue A and 
Kmin- For each A, Kmin is given by the largest zero of the function 2/2 (u)- 

Also in standard notations, the function 2/2 (m) can be written (apart from a 
constant factor) 

2/2(w) = u^+'^U Q - + i'^' 1 + «2i^; 
where U has the following definition: 



C/(a, b; u) 



sin nb 



M{a,b;u) u^~''M{a + 1 - b,2 - b;u) 



{a-by.{b-iy. 



(a -!)!(! -5)! 



The solution 2/2 (u), having a zero at m = Umin, is equivalent to requiring 

Arg M{u„itn) = -nn, 



(52) 



with n an integer. For very small Umm it is sufficient to include the leading 
term in an expansion of M . In this case eq. ( |5^ ) can be approximated by the 
relation 

mr 

nl Umin — ; 

V 

which means that n = 1 is relevant for the largest zero. A less severe approxima- 
tion (for finite Umm) is given by using the first three terms in the expansion of 
M(a, 6; u) and evaluating Arg[M(umm)] numerically. This leads to an approxi- 
mate relation between Kmin and the largest eigenvalue A, which agrees well with 
the numerical solution discussed in appendix |^. 
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C Random walk analogy 



As mentioned in the main text, for a running coupling constant the chain has 
similarities with a random walk in a force field. Let us first study a simple 
example to show what happens if a stochastic process of the Brownian motion 
type contains a "force". Consider a set of discrete "space-locations" indexed 
1, . . . , j, . . . and assume that at the discrete "times" t — 1, . . . , a, . . . there is a 
number of "objects", n(t = a,j) at the location indexed j. Assume further that 
there is a rule so that the number of objects changes with time according to 

n{a + l,j) = i-t-^n(a, j + 1) + ^-y^n(a,j - 1) (53) 

Thus there is for a positive a a tendency for the sites with a larger space index 
to lose to the sites with a smaller index as the time goes by, i.e. the change in 
the distribution contain a "force" directed towards small space indices. 

The equation can (by subtraction of n{a,j) on both sides and a limiting 
procedure in space x and time t) be rewritten into (the space-size and time-size 
parameters are d and r, respectively) 

dn d? d'^n ad dn 

For a = we obtain the well-known diffusion equation for n, i.e. a Gaussian 
with an average and a width both increasing with time ~ \/t. For a > we find, 
however, that starting with a narrow distribution at < = 0, the last term in Eq 
( p4[ ) becomes increasingly important when t increases. Thus the time derivative 
becomes reduced. For t > r/a^ the two terms on the right hand side become 
of equal magnitude in the dominating interval around x = d^JTfr. This implies 
that the solution approaches the time-independent distribution 

n ~ exp(— ax/c?) (55) 

This is just what happens for e.g. the atmospheric density of the earth. 

This example is not totally a relevant one for the BFKL mechanism because 
the occurring weight distribution has no simple probabilistic interpretation (al- 
though positive definite it is non-normalized as the weight corresponds to a den- 
sity). Nonetheless there are obvious similarities, all weights contain the notion 
of a "direction" , in this case a preferred direction towards small k- values in the 
cascade chains because of the running coupling, and the result is asymptotically 
an exponential distribution in k. 
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